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A  STREAMBED-PARTICLE  MODEL-STUDY  FACILITY 
USING  HYDROXYETHYLCELLULOSE 
SOLUTIONS  AS  A  FLUID 

By  Neil  L.  Coleman  and  Wilbert  M.  Ellis1 

ABSTRACT 

Variable-viscosity  long-chain  polymer  solutions  will  be  used  to  provide 
a  wide  range  of  Reynolds  numbers  in  a  new  facility  for  experimenting 
with  the  hydrodynamic  forces  exerted  on  streambed  particles.  Hydrox- 
yethylcellulose  solutions  tested  in  the  facility  proved  to  have  a  desirable, 
nearly  Newtonian  behavior  and  sufficient  viscosity  stability.  The  facility 
can  provide  the  needed  range  of  Reynolds  numbers,  while  maintaining 
flow-velocity  profiles  typical  of  streambed  conditions.  KEYWORDS:  hy- 
draulics, hydraulic  laboratories,  hydraulic  modeling,  sediment,  sediment 
transport,  soil  erosion. 


INTRODUCTION 

The  mechanics  of  the  beginning  of  streambed- 
particle  motion  are  important  in  problems  of 
stable  channel  and  riprap  designs.  Many  ex- 
perimental studies  have  been  made  to  investi- 
gate the  action  of  flow-induced  forces  on  a 
streambed  particle.  In  general,  where  the  parti- 
cle is  far  from  the  water  surface,  the  action 
of  streamflow  upon  an  alluvial  particle  is  a 
Reynolds-similitude  phenomenon.  For  this  rea- 
son, experiments  have  been  made  using  scaled- 
up  Reynolds  models  of  a  streambed  particle. 
In  these  simulations  the  character  of  the  local 
flow  and  the  application  of  forces  to  the  particle 
can  be  observed  in  a  detail  not  possible  in  flume 
experiments. 

In  model  studies  of  particles  in  the  silt,  sand, 
or  gravel  size  ranges,  obtainment  of  dynamic 
similitude  requires  operation  of  the  model  at 
low-particle  Reynolds  numbers.  Since  the  model 
particle  itself  is  large,  even  though  the  proto- 
type is  small,  and  since  the  velocities,  forces, 
and  so  forth  in  the  model  must  be  large  enough 


1  Geologist  and  hydraulic  engineering  technician, 
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search Service,  Oxford,  Miss.  38655. 


to  be  measured  conveniently,  low  Reynolds  num- 
bers must  be  obtained  by  means  of  working 
fluids  of  a  high  viscosity.  Furthermore,  since 
a  wide  range  of  Reynolds  numbers  is  generally 
required  in  the  model,  it  is  advantageous  to  be 
able  to  vary  selectively  the  viscosity  of  the 
working  fluid.  Early  work  with  a  large-scale 
particle  Reynolds  model  indicated  that  such 
conditions  might  be  accomplished  with  aqueous 
solutions  of  a  long-chain  polymer  as  a  working 
fluid,  if  a  polymer  with-sufficiently  Newtonian 
viscosity  characteristics  could  be  found  (3.) 2 

Many  long-chain  polymers  are  cellulose  com- 
pounds, several  of  which  are  used  industrially 
as  thickeners,  sizers,  and  adhesives.  They  are 
commercially  produced  in  a  variety  of  types. 
Depending  on  type,  they  can  be  dissolved  in 
water  to  produce  non-Newtonian  fluids  with 
drag-reducing  properties  or  nearly  Newtonian 
fluids  that  have  almost  stable  viscosities  over 
a  wide  range  of  shear  rates. 

Because  of  their  importance  in  naval  engi- 
neering and  other  applications,  the  drag-reduc- 
ing polymer  solutions  have  been  intensively 
studied  for  nearly  30  years  (6).  Perhaps  be- 

2  Italic  numbers  in  parentheses  refer  to  items  in  "Lit- 
erature Cited"  at  the  end  of  this  publication. 
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cause  of  the  emphasis  that  has  been  placed  on 
drag  reduction,  polymer  solutions  with  nearly 
Newtonian  characteristics  have  not  been  used 
as  working  fluids  in  model  studies  where  a 
wide  range  of  Reynolds  numbers  is  desired. 

This  paper  describes  a  streambed-particle 
Reynolds-model  facility  designed  to  operate 
with  Newtonian  polymer  solutions  to  study  lift 
and  drag  forces  on  streambed  particles,  flow 
configuration  and  pressure  distribution  around 
particles,  and  the  beginning  of  particle  motion. 
We  describe  the  construction  and  instrumenta- 
tion of  the  facility,  and  the  range  of  experi- 
mental flow  conditions  attainable  by  using  a 
nearly  Newtonian  polymer  solution. 

THE  MODEL  FACILITY 

Figure  1  is  a  schematic  of  the  particle  model- 
study  facility  that  has  been  built,  a  tunnel  that 
recirculates  water.  It  consists  of  (a)  a  rectan- 
gular upper  conduit  within  the  working  section 
of  the  model,  (b)  a  polymer  mixing  box,  (c)  a 
return  pipe,  and  (d)  a  propeller  pump  driven 
by  a  variable-speed  device  by  a  5-hp  motor. 
The  maximum  discharge  of  the  pump  is  1,100 
gal/min. 

The  water  tunnel  has  an  auxiliary  pump  and 
gravity  return  system  (e)  for  removing  fluid 
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Figure  1. — Schematic  of  the  model  facility:  (a)  rec- 
tangular conduit,  (b)  polymer  mixing  box, 
(c)  return  pipe,  (d)  propeller  pump,  (e)  auxiliary 
pump,  (f)  holding  tank,  (g)  drain,  (h)  air-vent 
valve,  (i)  water- filling  valve,  (j)  bleed  valve, 
(k)  circular  opening,  (m)  packing  gland, 
(n).  thermometer. 


and  storing  it  in  a  holding  tank  (f).  This  sys- 
tem permits  temporary  drainage  to  allow 
access,  through  appropriate  hatches,  to  the 
working  section  without  wasting  the  polymer 
solution.  For  changing  polymer  solutions,  there 
is  an  external  drain  (g).  The  polymer  mixing 
box  has  a  removable  lid  fitted  with  an  air-vent 
valve  (h)  and  a  water-filling  valve  (i)  with  a 
hose  fitting.  A  small  bleed  valve  (j)  in  the 
upper  conduit  upstream  of  the  working  section 
allows  withdrawing  fluid  samples  to  determine 
viscosity  and  density. 

To  use  polymer  solutions  economically  in  the 
model  facility,  an  initial  solution  of  a  relatively 
high  viscosity  must  be  mixed  and  a  complete 
series  of  experiments  performed  with  this  fluid. 
Then,  a  new  fluid  with  a  low  viscosity  is  pro- 

Table  1. — Properties  of  HEC  solutions  at  77°  F 


Concentration 
of  HEC 

(%) 


Density 
(slugs/ft3) 


0.2 


.6 


1.934 


1.934 


1.934 


1.934 


1.0 


2.0 


1.934 


1.953 


Shear  rate1 

Kinematic 
viscosity!  - 
(ft2/s) 

1  Qfi  v  1  Hi 

9  CIS  v  1  fi~5 

o.oo  X  1U 

9  A9  \y  1  H-5 

1        v  1  (15 

l.oo  X  J-^ 

1   M  v1(l-J 
l.oo  X  *U 

a  qq  yino 

O.OO  /\  IU 

4  ^vfi  V  1  fl-5 
*±.U\_>  X 

1.38  X101 

4.35X10"5 

3.42x101 

4.06  XlO"5 

6.88  xlO1 

4.19X10"5 

8.00  XlO4 

3.18X10-5 

3.44x10° 

1.16  XlO"4 

6.88x10° 

1.13  XlO"4 

1.38  xlO1 

1.10  XlO"4 

3.44  XlO1 

1.04  XlO"4 

3.54  x10s 

7.20x10-5 

6.90X10-1 

2.43  XlO"4 

1.72x10° 

2.22  XlO"4 

3.44x10° 

1.98  XlO"4 

6.88x10° 

1.94  XlO"4 

1.38X101 

1.95  XlO"4 

1.96x10* 

1.30  XlO"4 

6.90x10"! 

3.90  XlO"4 

1.72x10° 

3.79  XlO"4 

3.44x10° 

3.49  XlO"4 

6.88x10° 

3.46  XlO"4 

1.38x101 

3.48  XlO"4 

1.09  XlO4 

2.34  Xl0~4 

3.40x10-1 

2.99X10-3 

6.90x10"! 

3.08X10-3 

1.72x10° 

3.12x10-3 

1.61x10s 

2.18x10-3 

1  Last  value  in  each  set  of  measurements  was  de- 
termined with  a  falling-ball  viscosimeter,  other  values 
with  a  rotating-spindle  viscosimeter. 

2  Each  value  is  the  average  of  3  determinations  made 
at  the  same  shear  rate. 
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duced  by  controlled  dilution  of  the  original 
fluid,  followed  by  another  series  of  experiments. 
The  process  can  be  repeated  several  times,  so 
that  an  extremely  wide  range  of  experimental 
viscosities  can  be  obtained  with  a  relatively 
small  amount  of  polymer. 

The  upper  rectangular  conduit  (31  ft  long, 
0.916  ft  wide,  and  0.583  ft  high)  of  the  model 
facility  is  constructed  of  transparent  plastic. 
The  working  section  (a)  is  a  bed  of  i/o-in- 
diameter  plastic  spheres  glued  in  a  close-packed 
configuration  to  the  floor  of  the  conduit,  simu- 
lating in  scaled-up  form  a  small  streambed 
segment.  The  leading  edge  of  the  bed  of  spheres 
is  22.5  ft  from  the  upstream  end  of  the  conduit. 
The  bed  is  4.6  ft  long.  At  a  point  on  the  center- 
line  of  the  bed,  and  3.4  ft  below  the  leading 
edge,  a  3-in-diameter  circular  opening  (k)  in 
the  conduit  floor  allows  for  the  installation  of 
equipment  for  measuring  forces  on  a  bed  parti- 
cle. A  packing  gland  (m)  in  the  top  of  the 
conduit,  3.35  ft  from  the  leading  edge  of  the 
bed,  allows  a  total-head  tube  to  be  inserted  for 
determining  the  velocity  distribution  over  the 
bed.  For  monitoring  hydrostatic  pressure,  two 
piezometer  taps  (not  shown)  are  installed  in 
the  sidewall  of  the  conduit  near  the  openings 
(k)  and  (m).  A  thermometer  (n)  has  been 
installed  through  the  top  of  the  conduit,  just 
below  the  trailing  edge  of  the  bed  to  allow 
fluid  temperature  to  be  observed  during  an 
experiment. 


PROPERTIES  OF 
POLYMER  SOLUTIONS 

The  polymer  compound  selected  for  use  is 
hydroxyethylcellulose  (HEC),  which  is  com- 
mercially available.  According  to  the  manu- 
facturer, the  compound  is  nearly  Newtonian 
over  a  wide  range  of  shear  rates  and  relatively 
unsusceptible  to  viscosity  degradation  by  stir- 
ring or  pumping,  and  it  can  be  completely  dis- 
persed in  water  by  stirring  before  it  begins  to 
dissolve.  Such  characteristics  insure  formula- 
tion of  an  homogeneous  solution,  without  resid- 
ual clumps  of  partially  dissolved  polymer. 

Small  samples  of  HEC  solutions  at  various 
concentrations  were  analyzed  for  differences 
in  density  and  viscosity  (under  varying  shear 
rates),  to  determine  the  extent  to  which  these 


solutions  might  exhibit  non-Newtonian  be- 
havior. These  analyses  were  carried  out  at  a 
standard  temperature  of  77°  F  (25°  C).  The 
results  are  given  in  table  1.  Each  viscosity  value 
is  the  average  of  three  determinations  made 
at  the  same  shear  rate. 

The  solution  densities  (table  1)  were  found 
by  weighing  measured  volumes  of  the  different 
solutions.  The  table  shows  that  the  densities 
of  HEC  solutions  at  77°  F  do  not  differ  ap- 
preciably from  the  density  of  water  until  a  2  % 
concentration  is  reached.  The  kinematic  viscosi- 
ties in  the  table  were  measured  with  a  com- 
mercial rotating-spindle  viscosimeter  capable  of 
exerting  eight  different  shear  rates  on  a  sample, 
and  also  with  a  commercial  falling-ball  viscosi- 
meter. In  the  rotating-spindle  viscosimeter,  a 
cylindrical  spindle  with  radius  r  is  inserted  into 
a  small  sample  container,  so  that  the  fluid  is 
confined  to  a  restricted  annular  space  with 
radius  Ar,  small  relative  to  r.  The  shear  rate  (??) 
for  the  instrument  can  then  be  approximated 
as  a  function  of  spindle  speed  (<*>)  as 


where  Ut  is  the  tangential  velocity. 

With  the  falling-ball  viscosimeter,  the  ball 
falls  along  the  centerline  of  a  tube,  so  that  the 
fluid  sample  is  forced  through  a  small  annulus 
at  the  diameter  of  the  ball.  The  maximum  shear 
rate  (rj)  in  the  annulus  is  at  the  ball  surface 
and  can  be  approximated  from  considerations 
of  the  velocity  distribution  along  any  radius  of 
the  annulus,  and  shown  as  a  function  of  the 
ball  Reynolds  number, 


where  Rb  is  the  ball  radius,  W  the  ball  velocity, 
and  v  the  kinematic  viscosity  of  the  fluid. 
Apparently  (table  1),  the  falling-ball  viscosi- 
meter exerts  substantially  higher  shear  rates 
than  does  the  rotating-spindle  viscosimeter. 

The  non-Newtonian  behavior  of  the  HEC 
solutions  is  not  pronounced:  kinematic  viscosity 
for  the  six  solutions  in  table  1  decreased  about 
34%  over  an  average  100%  increase  in  shear 
rate.  The  variation  of  kinematic  viscosity  with 
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shear  rate  for  these  polymer  solutions  can  be 
represented  by 

HI 

where  vj  is  the  kinematic  viscosity  at  the  unit 
shear  rate,  i^.  For  a  perfectly  Newtonian  fluid, 
the  exponent  N  would  have  the  value  zero.  For 
the  polymer  solutions,  N  varies  with  concentra- 
tion, as  shown  in  figure  2.  The  relative-viscosity 
function  expressed  by  equation  3  is  depicted  in 
figure  3. 

Evidently,  there  are  two  ways  of  representing 
the  apparent  viscosity  of  a  polymer  solution 
used  as  a  working  fluid  in  the  model  facility. 
The  first  is  to  use  the  unit  shear  rate  (r^)  for 
a  particular  solution,  at  a  particular  fluid  tem- 
perature. However,  this  is  not  realistic  with 
regard  to  actual  shear  rates  expected  in  the 
model.  The  second  method  is  to  use  the  nearly 
constant  (quasi-Newtonian)  viscosity  exhibited 
by  the  polymer  solutions  at  high  shear  rates. 
We  adopted  this  method  to  study  the  flow 
conditions  in  the  model  facility,  since  the  high 
shear-rate  viscosity  of  a  given  working  fluid 
can  be  easily  obtained  with  the  falling-ball 
viscosimeter. 

PERFORMANCE  OF  THE 
FACILITY 

We  tested  the  series  of  HEC  solutions  for 
which  temperature-viscosity  functions  are  given 
in  figure  4,  to  verify  that  the  polymer  solutions 
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Figure  2. — Variation  of  the  exponent  N  with  polymer 
solution  concentration. 


would  maintain  constant  viscosity  during 
pumping  and  to  define  the  velocity  distributions 
and  the  range  of  experimental  conditions  that 
could  be  generated  in  the  facility. 

The  viscous  stability  of  the  HEC  solutions 
was  verified  by  comparing  the  temperature- 
viscosity  function  immediately  after  mixing 
and  after  several  hours  of  experimental  pump- 
ing. Figure  5  shows  the  comparison  for  a  0.3  % 
HEC  solution. 

Because  of  the  wide  range  of  viscosity  of  the 
HEC  solutions,  the  flows  generated  in  the 
facility  can  be  expected  to  display  time-mean- 
velocity  distributions  ranging  from  the  para- 
bolic distribution  of  a  completely  viscous  flow 
to  the  logarithmic  profile  of  a  turbulent,  rough 
conduit  flow.  These  velocity  distributions  can 
all  be  accommodated  (7)  on  a  single  graph 
plotted  in  the  coordinates 

^!H(¥)+T»(l).  «> 

where  U  is  the  time-mean  velocity,  y  the  dis- 
tance from  the  velocity  measurement  to  the 
virtual  origin  of  the  profile,  £7*  the  shear  ve- 
locity at  the  virtual  origin,  «  the  Karman  con- 
stant, and  8  the  boundary-layer  thickness,  taken 
as  the  distance  from  the  point  where  the  maxi- 
mum velocity  Um  occurs  in  the  conduit  to  the 
virtual  origin.  The  ratio  a  £7/ £7*  is  the  velocity 
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Figure  3. — Relative-viscosity  function  for  HEC  solu- 
tions. 
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estimated,  using  an  adaptation  of  the  cross-plot 
method  (8). 

In  figure  6,  the  velocity  shift  AU/U*  is  plotted 


2.0 
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Figure  4.— Viscosities  of  working  fluids  for  solutions  of  various  nominal  concentrations 

shift  caused  by  boundary  roughness,  and  the 
function  (II/k)  l»(y/8)  ]  is  Coles'  wake  func- 
tion U),  in  which  n  is  the  coefficient  of  wake 
strength.  The  function  f(U*y/v)  is  the  "stand- 
ard" velocity-profile  function  that  can  be  de- 
rived from  simple  dimensional  analysis  or  from 
mixing-length  theory  (8). 

Velocity  profiles  were  measured  in  the  model 
facility  with  an  appropriately  designed  and 
calibrated  total-head  tube  (2)  and  a  differ- 
ential pressure  transducer.  The  total-head  tube 
was  inserted  through  the  packing  gland  (fig.  1, 
m).  It  was  constructed  so  that  its  tip  would 
traverse  a  vertical  located  3.2  ft  from  the  lead- 
ing edge  of  the  test  bed,  or  just  ahead  of  the 
transducer  insertion  port  (fig.  1,  k).  In  making 
the  velocity  measurements,  we  took  the  virtual 
origin  of  the  profile  as  one-fifth  of  a  sphere 
diameter  (0.2Z)  or  D/5)  below  the  top  plane 
of  the  bed,  as  recommended  by  Einstein  and  T'°F 
El-Samni  (5). 

The  velocity  profiles  were  plotted  in  U/Um       FlGV™  5— Viscosity-temperature   curve   for   a  0.3% 

(Umy/V)  coordinates,  and  U*  and  At// U-  were  !?Uti°n  -(a)  after  initial  mixing'  (b)  after 

,  a,nu  au/u.  wwc  7  hours  of  pumping  during  experiments. 
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-4.2, 


(6) 


which  is  identical  to  that  for  sand  roughness, 
except  the  constant  is  4.2  rather  than  3.6,  indi- 
cating that  a  bed  of  closely  packed  spheres  is 
slightly  smoother  than  a  bed  of  natural  sand 
particles.  The  above  results  agree  generally 
with  those  of  other  experimenters  (10). 

The  velocity  profiles  are  plotted  in  figure  7. 
The  functions  included  in  the  figure  are,  except 
for  the  introduction  of  A~U/U#,  the  usual 


U  _  U*y 
V.,  V 


(7) 


Figure  6. — Velocity-shift  function, 
against  the  bed-shear  Reynolds  number: 


for  the  viscous  sublayer,  and  the  Karman- 
Prandtl  law  of  the  wall, 


R, 


U,D 


(5) 


U 

IT, 


2.303 


-  log 


+  4.9, 


(8) 


Figure  6  shows  that  bed  roughness  begins  to 
affect  the  flow  at  U«D/v  equal  to  approximately 
4.  The  transition  regime  ends,  and  the  fully 
rough  flow  begins  at  about  50.  In  the  fully 
rough  regime,  the  function  is 


for  the  inertial  region.  For  plotting  figure  7, 
the  Karman  constant,  «,  in  equation  8  was 
taken  as  0.41.  The  dash  curve  in  the  figure  is 
the  usual  empirical  transition  function  between 
the  viscous  and  inertial  regions.  The  velocity 
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Figure  7. — Typical  velocity  profiles. 


6 


profiles  labeled  A  in  the  figure  are  characteris- 
tic of  totally  viscous  flow  and  have  the  equation 

which  is  expanded  to 

for  plotting  in  figure  7.  These  profiles  are 
asymptotic  to  the  viscous-sublayer  function. 
The  profiles  labeled  B  are  typical  of  flows 
where  the  displacement  Reynolds  number, 


(ID 


is  less  than  725  (8).  These  flows  have  a  strong 
wake  component  and  no  logarithmic  region. 
Profiles  labeled  C  are  typical  of  flows  with 
Rg*  greater  than  725,  with  significant  intensity 
of  turbulence  in  the  inner  part  of  the  conduit. 
These  profiles  show  both  a  logarithmic  part  and 
a  wake  region.  As  conduit  turbulence  increases, 
the  wake  component  decreases  (4),  and  pro- 
files like  D  appear.  These  profiles  have  a  loga- 
rithmic region  but  no  wake  region. 


The  bed-shear  Reynolds  number  (R*)  is  the 
independent  variable  on  which  the  velocity  re- 
duction in  figure  6  depends.  This  Reynolds 
number  is  an  indicator  as  to  whether  the  flow 
over  the  bed  is  hydraulically  smooth  or  rough, 
and  whether  or  not  a  particle  on  the  bed  is 
experiencing  viscous  or  turbulent  approach 
flow.  Thus,  R*  is  one  type  of  particle  Reynolds 
number,  a  second  is  the  conventional  form 


where  U0  is  the  approach-flow-velocity  vector 
acting  along  a  line  through  the  particle  center. 
This  Reynolds  number  characterizes  the  local 
flow  around  an  isolated  streambed  particle  and 
serves  as  the  primary  independent  similitude 
parameter  in  lift  and  drag  coefficient  studies. 

The  joint  variation  of  R;  and  Re  is  plotted  in 
figure  8.  Also  included  for  comparison  are 
the  ranges  of  Rr  obtained  by  Coleman  (3)  and 
Aksoy  (1),  and  the  ranges  of  R*  obtained  by 
Watters  and  Rao  (9).  Apparently,  the  facility 
described  here  has  broader  experimental  po- 
tential than  those  used  in  the  other  model 
studies.  The  capability  of  operating  over  a 
range  from  totally  viscous  flows  to  fully  tur- 
bulent rough  boundary  flows  is  particularly 


Figure  8. — Range  of  conditions  attainable  in  the  model  facility. 
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advantageous  for  streambed-particle  model 
studies. 


By  using  a  water-tunnel  facility  and  HEC 
solutions  with  a  wide  range  of  viscosities  as 
working  fluids,  researchers  can  conduct  stream- 
bed-particle  model  studies  over  a  wider  range 
of  Reynolds  numbers  than  obtained  previously. 
Particle  Reynolds  numbers  from  1.0  to  10,000 
can  easily  be  obtained,  and  flow  conditions  can 
range  from  completely  viscous  to  completely 
hydraulically  rough. 
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